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1 Introduction 


There are two main lines of research of chaotic behavior in General 
Relativity: one deals with chaoticity associated to inhomogeneous cosmo¬ 
logical models, as in Bianchi IX model [|l|. The invariance of the relativistic 
theory challenges the proper concept of chaos in this approach Although, 
a gauge invariant method based on Maupertuis principle is proposed in [Q 
and applied with little changes to the Bianchi IX model |^. The method 
points to the occurrence of local instability in this cosmological model near 
its singularity. This adds to a plethora of early analysis and numerics (see 
for instance the references cited in [Q), which reinforces the increasingly ac¬ 
cepted idea that Bianchi IX is chaotic in some meaningful sense. 

The other line assumes a given geometry and looks for chaotic behav¬ 
ior of geodesic motion in this background. In this case, the geometry can 
be taken as either an approximate or an exact solution of Einstein equa¬ 
tions. Examples of chaotic geodesic motion are considered in p - P] (exact 
geometry) and in p, |I^ (approximate). 

Due to its universality and intrinsic mathematical interest, models in 
which unstable periodic orbits (UPOs) are subjected to small periodic per¬ 
turbations has been one of the main paradigms of deterministic chaos [|^ . 
An analytical tool to study these models is the Melnikov function that de¬ 
scribes the transversal distance between the unstable and stable manifolds 
emanating from an UPO. Its isolated odd zeros indicate the crossing of these 
manifolds, hence the onset of chaos [^, |T^. Examples of applications of the 
Melnikov method in gravitation are: the gravitational collapse of cosmologi¬ 
cal models |^, the study of orbits around a black hole perturbed by either 
gravitational radiation or an external quadrupolar shell []^. Also for par- 
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tides moving in several models of attractive centers periodically perturbed 
this method has been applied 0- 

In this work we study the occurrence of chaos for test particles moving 
around a Schwarzschild black hole perturbed by a special class of gravita¬ 
tional waves 1^. Since this perturbation is an exact solution of the Regge- 
Wheeler equations (see for instance 0) we can study the range of the pa¬ 
rameters of the model in order to have chaos. In a similar situation is 
studied but the perturbations, albeit more general that ours, are known only 
in the high frequency limit. 

In the next section we present a summary of the Melnikov method. In 
Sec. 3 we review the homoclinic orbits for Schwarzschild solution, and present 
an explicit solution t = t{r) for this orbit. In sections 4 and 5 we consider the 
perturbations and the the application of the Melnikov method, respectively. 
In the last section we discuss some of the previous results and also compare 
our work with the application of the Melnikov method presented in P]. 


2 The Melnikov method 


Let us consider a Hamiltonian, 


p 




( 1 ) 


that admits at least one UPO with the corresponding homoclinic orbit, and 
also a small periodic perturbation described by the Hamiltonian function 
eHi{q,p,t). Then the transversal distance, in phase space, between the per¬ 
turbed unstable and the perturbed stable manifolds emanating from the UPO 


is proportional to [13 




( 2 ) 
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where the integral is taken along the unperturbed homoclinic orbit, to is an 
arbitrary initial time, and 


if T 

dpdq dqdp' 


( 3 ) 


with {q,p) being canonical conjugated variables. If there is an isolated in¬ 
tersection for some to, i.e., an isolated odd zero of M(to), then there will be 
one for every to- This inhnitely many crossings of manifolds will produce a 
tangle that is the signature of the homoclinic chaos [[ 

In particular, if Hi has the form 


la 


Hi{p, q, t) = Hi{p, q) cos(a;t). 


( 4 ) 


we have 13 


M(to) = cos(a;to) / {i/o,cos(a;t)(it-|-sin(a;to) / {Hq, Hi} sm{ut)dt 

J —oo J —oo 

( 5 ) 

Then, in the generic case, we will have isolated zeros when at least one of 
the integrals of the previous formula is different from zero. 


3 The homoclinic orbits 

We shall consider a relativistic particle moving in a £x spacetime described 
by the metric Qab- The world-line of a particle will be denoted by x°'{s) and 
its mass by fi. The motion of the particle can be obtained from the action 

S[x] = ^ J gahX°'x^ds. ( 6 ) 

This action is not reparametrization invariant, but is the simpler for our 
purposes, s is the proper time along the world-line. The canonical conju¬ 
gate momentum to is Pa = and satishes the mass shell constraint 
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g°-^PaPb = —/i^. The Hamiltonian of the system is 

H = (7) 

The background metric will be considered as the metric of a non rotating 
black hole, i.e., the Schwarzschild metric, 


ds'^ = —fdT^ + f ^dR^ + R?{d§‘^ + sin^ iddip'^), (8) 

where / = 1 — 2M/R. Since, T and ip are cyclic variables we have the 
conserved quantities: 

P _ fdt 

E = -pT = 

ds 

L = p^ = pR^ (9) 


For planar motion 'd = 7r/2 {p^ = = 0) we have the equivalent one 

dimensional problem of a particle of mass 2p^ moving in the the potential V 
and constant energy E'^ — p‘^, 


2Mp^ , 2ML^ 

It is convenient to work in dimensionless variables. Dehning r 
U = {2M/LfV we hnd 


( 10 ) 
R/2M and 


U 


( 1 -/ 3 ^) I 1 1 

/^3 


( 11 ) 


where (3 = — 12M‘^p‘^/L‘^ < 1. Since U{oo) ~ —0(l/r), in order to 

have a local maximum (or unstable point, r = r„) with I7(r„) < 0, we need 

0<f3< 1/2. 
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In Fig.l we plot U as a function of r = R/2M for various values of jS. 
The unstable periodic orbit (UPO) correspond in phase space to the point 
{r,pr) = (r„, 0). Now by taking f/(r„) = (2M/L)^(E^—/i^) we get the motion 
equation for the homoclinic orbit, 

dr _ _^ wp{r - l)(r - r„)(r^ - r) 
dt r3/2 ’ 

where = 3/(1 — 2/3), and = | y ^ 3(1 — 2/3/(2 — (5Y. We have used 
the constant of motion E to change the parameter of evolution s to the 
dimensionless time coordinate t = T/2M. So the homoclinic orbit is limited 
Iqj r = Tu and the maximum value r = (turning point). This motion 
equation admits the primitive. 


±W0t = [r{rm - r)Y/‘^/rm + (2 + 2r„ + r^) arctan 


2 ^ 5/3 

rm/r - 1 

- ru)y/rm- Tu > 

rm/ru - 1 

2 ^ / 

Tm-r 


r{rm - 1 )’ 


r)/r 


( 13 ) 


We have chosen the constant of integration to have t{r = r^) = 0. In Fig.2 
we show the positive branch of ([T^) for different values of /3. The particle 
takes an inhnite amount of time to depart (arrive ) from (to) the homoclinic 
point where the UPO is located, this is a universal behaviour for these type 
of orbits. The explicit expression (|^) for the homoclinic orbit will play an 
important role in our study of the Melnikov function and its use will be a 
signihcant departure from the treatment given in wherein other interesting 
graphics related with the homoclinic orbit are presented. 
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4 The perturbations 


The perturbations of the black hole that we shall consider are of the particular 
class gab + ehab with gab given by (H) and 

hTT = -fXPi cos(aT) 
hnR = f-^YPi cos{aT) 

= R^{ZPi + WPi,^ ) cos(crT) 
hip^p = P^ sin^ '&{ZPi + WPi,^ cot "d) cos(crT), (14) 

where X, Y, Z and W are functions of R to be determined by the Einstein 
equations Rab{gcd + ^hcd) = 0. Pi = Pi(cos'd) are the usual Legendre poly¬ 
nomials. When the above mentioned Einstein equations are expanded up to 
the first order in e we find a linear system of first order differential equa¬ 
tions for the functions X, Y, Z and W [j^. The particular separation of the 
angular part of the perturbation (|T^) was found by Friedman |]^ . Combina¬ 
tions of this variables reduce these equations to the Zerilli equation, i.e., to 
a Schrodinger type of equation with a nontrivial potential . A particular 
solution of the hrst order differential equations is 


X = pq, Y = 3Mq, Z = {R — 3M)q, W = Rq, 


( 15 ) 


where 


q = f^l'^lR\ p = M- 


+ a^R^ 


( 16 ) 


R-2M 

This solution was found by Xanthopoulos [|T^] using an ad hoc method. We 
note that it can be easily found in the Zerilli formulation and corresponds to 
zero Zerilli function ZY) ^ 


= 


nR + 3M R 


(—W-E), 


( 17 ) 
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where n = {I - l){l + 2)/2. 

From Eq.([^) we have for the perturbed system 

- eh-^aVb = -/i^ ( 18 ) 

For an even I we have that the perturbations ( 0 ) are even functions in 
the variable cos'd, i.e., we have perturbations with reflection symmetry with 
respect the plane '& = 7t/2. Then, for a particle moving in the plane '& = 71 /2 
these type of perturbations do not change the plane of the orbit. In other 
words, the perturbation introduces a “perpendicular force” that is an odd 
function of cos'd. We hnd for particles moving in the plane'd = 7r/2, 


-pt = Ho + ei/i. 


with 

E 1 f 

Hi = -[yhTT - ^{fpRhRR + 

and 


hrr = -fpqNi cos(aT), 
hRR = —^Wcos(crT), 

= R^{R — 3M)qNi cos(crT), 


where Ni is zero for odd I and 


Ni = 


1 ■ 3 ■ ■ ■ (2n - 1) 

2■4 ■ ■■2n 


(/ = 2n). 


( 19 ) 


( 20 ) 


( 21 ) 


( 22 ) 


The order of Legendre polynomial I enters in a trivial way in this class of 
planar perturbations, as a multiplication by the constant W- The frequency 
a not only appears in the oscillating function cos((jT) but also in p. 





5 The Melnikov function 


Now we shall use the fact that we know explicitly the function t = t{r) for 
the homoclinic orbit to evaluate the Melnikov function. We hnd 


M{to) = — sm{ujto)K, 


with 


K = ^NiJ, 

M 

m 

J= F{r, (3) sm{(jt{r))dr, 

J ru 

where uj = 2Ma is a dimensionless frequency and 
F{r,i3,uj) = 


1 I _ / , 27/^ 1 . Vy _ 

2^2 « 2 r 5 4 r^( 2 M )2 2r^ {2My 


'y h 

iJ 


+ 77w(l ~ ‘^'yf'^)hrr — i{~^hrr + —hr 


2 r 2 


dt 


hit — 


2 r 2 

3 


1 - 


1 + 4a;r^' 

2 r/ , 


2r2/V2 ’ 


K 




{2My 


= f/2(2r-3)/2. 


(23) 

(24) 

(25) 


(26) 


The constant 7 is dehned in terms of as 7 = ^(1 + (3){2 — /d)^, and 


(27) 


By t = t{r) in (^) we mean the positive branch of (1^). We have used the 
homoclinic orbit to map the inhnite interval of the integration to a hnite one. 
Also, all the quantities appearing in the dehnition of J are dimensionless. 

Therefore, to have the homoclinic tangle of orbits we need J 7 ^ 0. 
The exact computation of J even in the simple case of the perturbations 
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under consideration is hopeless. Nevertheless, the integrand of J is explicitly 
known in terms of the variable of integration. This fact can be used to study 
the zeros of J = J{P,uj), or better the range of f3 and u in which J 7 ^ 0, in 
a simple graphic way. 

Firstly, it is instructive to plot the function S = sin (cut (r)) for a 
hx f3 and different values of 00 . We hnd (see Fig.3) that, depending on 
the value of cn. S' does not change sign in a large portion of the interval 
f'uiP) ^ r < TmiP), e.g., the two bottom curves (from right to left). We also 
have a very rapid oscillation near the position of the UPO, that is not shown 
in the graphic. The values of the parameters are u = 0.15 (top, from right 
to left), 0 . 1 , and 0.05 (bottom) and /? = 1/4. Also it is illuminating to make 
a graphic of F{r, (3, u) for different values of (3 and uj in the range of interest 
< r„i{(3)- Fig. 4 shows that this function is very well behaved and 
does not change of sign in the interval of interest. The integrand of J, i.e., 
F{r, P,u) sm{ujt{r)) with f3 = 1/4 and a;=0.1 is presented in Fig.5. The area 
under the curve is clearly different from zero, then J 7 ^ 0. In this case the 
perturbation will originate a chaotic motion of the test particle. The value 
cu = 0.1 was chosen from the relation 

~ 7r/4, (28) 

where is the value of the homoclinic orbit that correspond to the maximum 
value of its curvature k = t"/{I + (F)^)^/^. For /? = 1/4 we have 4=29.3691 
and ci;fc=0.1069. Therefore, we can say that the perturbations with oj < uik 
will originate chaotic behaviour. 
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6 Discussion 


To better understand the particular perturbation (|I4|) we have computed 


the Kretschmann scalar for the perturbed metric for the particular values of 
I = 2 and Z = 3, we find 



(29) 


the numerical constant ki is ^2 = 58 and = —3/2. Thus, we have a 
Kretschmann scalar decaying as R~^ for R 3> 2M, indicating that part of 
the wave has origin in the black hole. 

In this work we have shown that for a particular class of gravitational 
perturbations we have that the particles moving along timelike geodesics will 
present a chaotic motion when the frequency of the perturbation is a; < 

In the complementary case, uj > ujk our simple graphic method does not 
apply. But, anyway one can compute numerically the integral J with a great 
precision ( the integrand is known explicitly). We found that for uj > ook, 
J has positive and negative values depending on the values of (3. Then 
for certain specihc values of (3 we will have J{j3,uj) = 0 and the Melnikov 
method does not apply to these cases. From the previous analysis it does 
look like that the generic situation is chaotic, i.e., only for a numerable set 
of frequencies a;, given a value of (3, we should not have chaos. The number 
of elements of this set increases with a;, but is still numerable. To be more 
precise, the Melnikov method fails to predict chaos only for a a numerable 
set of values of the parameters. We note that seldom one can find so easily 
the parameter range of a given problem that will produce a chaotic situation. 
Moreover, for parameters outside this range it is reasonable to guess that only 
for a numerable set of frequencies the integrability is preserved. Therefore, 
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based in our previous analysis, we can conjecture that chaos is generic for 
the class of perturbed systems studied. In our case, the limit lim^^^oo 
is not dehned. Note that in F{P,uj) appears o;^ and in the integrand we also 
have cos cut. In the case analyzed in p the perturbation are quite general, 
but they are known explicitly only in the limit a; —cx). Their integral 
equivalent to our J has ±cx) integration limits and in the integrand they have, 
essentially, an oscillating function with frequency uj multiplied by another 
function proportional to uj. Since, in 0, every thing is known only in the 
limit of inhnite frequency, in order to predict chaotic behaviour, it is necessary 
to show that their integral is different from zero in the same limit. This point 
is unclear in the quoted paper. 

The authors thank CNPq and FAPESP for hnancial support and Prof. 
D.V. Gal’tsov (Moscow) for enlightening discussions about black hole per¬ 
turbations. 
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FIGURE CAPTIONS 


Fig.l. The effective U potential is plotted for the dimensionless parameter 
/3 = 1/3 (top cnrve), 1/4, and 1/5 (bottom cnrve); the maximnm value is 
located at = 3/(1 + /3). 


Fig.2. A graphic of the positive branch of Eq. (^) is pictured for values of 
the parameter /3 = 1/3 (top curve), 1/4, and 1/5 (bottom curve). We have 
that the particle takes a finite time to travel from to the vicinity of and 
then takes an infinite time to arrive to Tun itself, wherein is located the UPO. 


Fig.3. A graphic of 5 = sin[a;r(r)] for /? = 1/4, and oj =0.05 (bottom curve, 
from right to left), 0.1, and 0.15 (top curve). 


Fig.4. A plot of F{r,j3,uj) for (3=1/4:. 

Fig.5. A plot of F{r,/3,uj) sm{ut{r)) for (3=1/4, and uj= 0.1. We have a 
rapid oscillation near (=12/5 ) (not shown in the hgure), that produce 
very hne spines in the range [-0.03, 0.03]. 
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